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Propylene carbonate (PC) wets graphite with a contact angle of 31° at ambient conditions. Molecular 
dynamics simulations agree with this contact angle after 40% reduction of the strength of graphite-C atom 
Lennard- Jones interactions with the solvent, relative to the models used initially. A simulated nano-scale PC 
droplet on graphite displays a pronounced layering tendency and an Aztex pyramid structure for the droplet. 
Extrapolation of the computed tensions of PC liquid- vapor interface estimates the critical temperature of PC 
accurately to about 3%. PC molecules lie flat on the PC liquid- vapor surface, and tend to project the propyl 
carbon toward the vapor phase. For close PC neighbors in liquid PC, an important packing motif stacks 
carbonate planes with the outer oxygen of one molecule snuggled into the positively charged propyl end of 
another molecule so that neighboring molecule dipole moments are approximately antiparallel. The calculated 
thermal expansion coefficient and the dielectric constants for liquid PC agree well with experiment. The 
distribution of PC molecule binding energies is closely Gaussian. Evaluation of the density of the coexisting 
vapor then permits estimation of the packing contribution to the PC chemical potential, and that contribution 
is about 2/3rds of the magnitude of the contributions due to attractive interactions, with opposite sign. 



I. INTRODUCTION 

Unique properties of nanomaterials, and specifically of 
electrochemical double-layer capacitors (EDLC) based on 
carbon nanotube (CNT) forests,'^ arise from their large 
surface areas. Accurate descriptions of interfaces is an 
important challenge for modeling. Here we considered 
propylene carbonate (PC: FIG.jl]) as a solvent for electro- 
chemical double-layer capacitors, reporting experimental 
and molecular simulation results on PC interfaces. This 
validation is preparatory to direct simulation of dynam- 
ical filling and performance of CNT-based EDLCs. 

Differential filling of a charged nanotube forest from 
a bulk electrolyte solution is expected to be sensitive 
to the balance of attractive intermolecular interactions. 
A possible concern for the first simulations of CNT- 
based EDLCs was that models of the PC solvent had 
not been parameterized to describe contact with car- 
bon electrodes.!^ An initial molecular dynamics calcu- 
lation showed complete spreading of the modelled liquid 
PC on graphite,^ in disagreement with experimental ob- 
servation of a PC droplet on graphite reported below. 
That experience suggested the advantage of studying 
PC interfaces in validation of simulation models for elec- 
trochemical applications. The present work adjusts the 
Lennard- Jones interactions with the graphitic carbon to 
agree with the observed contact angle. 

Propylene carbonate is a non-reactive, low-toxicity, 
aprotic, highly polar dielectric solvent widely used 




FIG. 1. Propylene carbonate (PC) enantiomer with atom 
labelling used in this paper. 



in nonaqueous electrochemical systems. Simulation 
of interfaces of non-aqueous solvents for electrolytes, 
acetonitrile^ ^ and propionitrile,^ have bee n pre viously 
carried-out. Helpful molecular simulationP^^^ of the 
uniform liquid propylene carbonate solvent are also avail- 
able. A variety of molecular dynamics results are dis- 
cussed below, and technical details of the methods are 



collected in the Appendix (Sec IV). 



II. RESULTS 
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A. Droplet-on-graphite contact angle 

The observed contact angle (FIG. [2| is acute. This 
indicates good wetting and favorable PCigraphite inter- 
actions, though not the complete spreading that was 
obtainecP from simulation with initial models Trial- 
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FIG. 2. PC droplet on graphite. Contact angle measurements 
of 5 /il propylene carbonate drops on a pyrolytic graphite 
disc (Grade: PG-SN, Graphitestore.com) were performed un- 
der ambient conditions using a Rame-Hart contact angle go- 
niometer. The graphite disc was sheared against a paper wipe 
(Kimwipe) to reveal a fresh surface. The surface was then 
rinsed with water and allowed to air dry. A static contact 
angle of 31.4zbl.6° was obtained from 3 trials. 
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FIG. 4. The density, relative to the coexisting liquid density, 
in the central plug of material indicated by the shaded region 
of FIG. [3] Atom layers nearest to the graphite are concordant. 
The density oscillations of CI and C4 atoms conflict at larger 
distances from the graphite, and C4 prevails furthest from the 
graphite. 
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FIG. 3. Droplet silhouette, suggesting an Aztex pyramid! 
The indicated contact angle is found by fitting a straight line 
to the region 3.0 nm < r < 4.2 nm. The shaded region 
identifies the cylinder used to investigate the layering of the 
mean density in FIG. [4] 



and-error (Sec. IV A) found a simulation contact angle 
of about 31° with a 40% lowering of the strength of 
Lennard- Jones interactions associated with the C atoms 
of the graphite substrate (FIG.jsj). 

Simulations of water on graphite surfaces show that 
changing the number of water molecules from 1000 to 
over 17,000 has only a small effect on the contact angle. 
A droplet size dependences of such results are often as- 
cribed to a line-tension effect. The fine structure ob- 
served here in the three-phase contact region here would 
make definition and further investigation of a triple line 
complicated. Together with the distinctly non-spherical 
molecular shape and packing pattern (studied below), 
interesting structural transitions of PC molecules near 
charged electrodes are expected; those issues should be 
the subject of future studies. 

The millimeter length scales of the experimental obser- 
vation differ by 10^ from nanometer length scales treated 
by the calculation. Pore radii for ECDL capacitors 
are in the nanometer range. Therefore, the nano-scale 
fine-structure of the molecular dynamics result (FIG. [4| 
should be relevant to the anticipated studies of EDCLs. 




FIG. 5. Slab geometry, upper configuration from T=300K 
simulation, lower configuration from T=600K simulation. 



B. Liquid-vapor interfacial tensions 



Surface tensions were evaluated from simulations in 
slab geometry (FIG.jsj) at several temperatures. The PC 
liquid- vapor interfacial tensions (FIG.[6| agree accurately 
with experimental values (TABLE [ij) at moderate tem- 

fr-T \^^^ 

peratures. The simple formula 7(T) ~ 70 ( ^ \ , 



utilizing an approximate estimate of the critical expo- 
nent, extrapolates to 7 ^ at ^ 740K (FIG. [6|. 
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TABLE I. Experimental surface tension at 20C is 41.1 dyne/cm. Vapor densities pvap were obtained from FIG. 13 utilizing 
WHAM calculations at the lowest temperature. 



T(K) 


piiq (nm~^) 


Pvap (nm ^) 


-/cTln (^^j (kcal/mol) 


PvapkT (bar) 


p (barj^ 


7 (dyne/cm) 


300 


7.11 


1.6e-6 


-9.2 


6.2e-5 


7.6e-5 


40.8 


400 


6.46 


4.8e-4 


-7.6 


2.6e-2 


2.8e-2 


29.9 


600 


4.93 


7.0e-2 


-5.1 


5.8 


4.9 


9.7 
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FIG. 6. Surface tensions 7 = 

I / {pzz{z) - \ \pxx{z) ^pyy{z)WAz. Tlic cxperimeutal 
critical temperaturj^ is Tc = 763 K, whereas the triple 
temperature is about Ttripie ~ 220 K. 



FIG. 7. The carbonate plane lies fiat on the surface, project- 
ing C4 from the liquid toward the vapor phase. 



C. Liquid-vapor interfacial structure 

For the lowest temperatures considered here, the in- 
terfacial profiles of the atomic densities are distinctive 
(FIG. [7|): non-monotonic except for the density of the 
propyl carbon (C4). The coincidence of the positive lip 
in atom densities other than C4 with the deficit in the 
C4 density suggests that the PC molecule lies approxi- 
mately fiat on this interface while projecting C4 further 
toward the vapor. Direct interrogation of the orienta- 
tions of molecules with zq\ >1.5 nm (FIG. [t]) confirms 
this view (FIG.|8|. 



D. Structure of the coexisting liquid 



2.0 




The atom-atom intermolecular radial distribution 
functions (FIGS. 9] and 10) supplement the view that 
an important packing motif stacks carbonate planes of 
close PC neighbors with the outer (01) oxygen of one 
molecule snuggled into the positively charged propyl end 
of another molecule so that neighboring molecule dipole 
moments are approximately antiparallel. Such an ar- 
rangement is similar to the known crystal structure of 
ethylene carbonate,^ and to the PC dimer structures 



FIG. 8. For molecules with zci >l-5 nm (FIG.[7|, the proba- 
bility density for projection of unit vector normal to the car- 
bonate (-CO2-) plane onto the z-axis, perpendicular to the 
interface. In this interfacial layer, the most probable orien- 
tation aligns the carbonate plane parallel with the plane of 
the interface, with the C4 methyl group extended toward the 
vapor phase. Uz > 0.5 {0 < 60°) for about 50% of interfacial 
PC molecules. 
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FIG. 10. Taking CI as the center of the PC molecule, this 
characterizes the packing of centers, with molecular coordi- 
nation numbers of 10 - 12. The solid blue curves are the 
radial distribution functions for the closest, and 2nd-closest 
CI neighbors of a CI atom, peaked at 0.38 nm and 0.44 nm, 
respectively. The dotted curve is the estimate of the radial 
distribution of the closest CI neighbor of a CI atom based 
on the Poisson approximation.^ The dashed blue curve is 
the probability, with median point of about 0.52 nm, that a 
central CI atom has more than two (2) CI neighbors within 
that radius. The closest Cl-Cl contacts are associated with 
stacking, with substantial disorder, of carbonate planes as 
suggested by the embedded molecular graphic. 



FIG. 9. Structured close-contacts involve negatively charged 
01 with the opposite end of the PC molecule. 



established by electronic structure calculations.^ 



On the basis of the ClCl joint distribution (FIG. 10), 
the nearest two neighbors of a CI are physically distinct, 
imperfectly stacked on top and bottom of the carbonate 
plane. The closest CI neighbor of a CI atom is most 
probably at a radial displacement of 0.38 nm (FIG. 10). 



The observed distribution of the closest CI neighbor of 
a CI atom (FIG. [lO| is more strongly peaked than the 
simple Poisson-based estimate.l^ 

The dielectric constant of uniform PC liquid implied 
by these simulations (FIG. 11) agrees satisfactorily with 
experiment. 
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E. Balance of excluded volume and attractive interactions 

Thermodynamic energies for liquid PC may be ana- 
lyzed on the basis of quasi- chemical theor}^^^^ in which 
the paramount goal is molecular-scale physical clarity 
from thermodynamic characteristics. The quasi-chemical 
approach focuses on binding energies, e = U {N) — U {N — 
1) — /7(1), of individual molecules and introduces a con- 



FIG. 11. Dielectric constant of model propylene carbonate, 
evaluated following standard simulation methods.^ The cal- 
culation treated 600 PC molecules under periodic boundary 
conditions. Resulting values were averaged from the 40 ns 
production trajectories at constant pressure of 1 atm. The 
solid line is the fit of Cote, et al}^ to their experimental re- 
sults, the solid dot their measured result at 25C. Numerical 
values are given in TABLE \n\ 
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ditioning based on definition of an indicator function,^ 
X, so that X = 1 defines a logical condition or constraint 
that physically simplifies that statistical thermodynamic 
problem. 

For example, if close neighbors set particularly impor- 
tant or complicated or chemical interactions, then the 
logical X = 1 can indicate the absence of neighbors in a 
specified local region. I ndeed, for numerous applications 
to liquid water ,l ^^^ l ^^tf^^' x = 1 indicates that there are 
no 0-atoms of bath molecules within a specified radius 
from the 0-atom of a distinguished water molecule. Then 
X = if any solvent molecule is present in that inner- 
shell. For the application to CF4(aq)^^, x = if ^iny 
water molecule is in defined van der Waals contact with 
the polyatomic CF4 solute. 

With such a indicator function x specified, the partial 
molar Gibbs free energy in excess over the ideal contri- 
bution, /i^®^^ = M ~ can be expressed as 

/3M(-) = -ln({x))o 

+ \nj e''^P(e|x = l)de 

+ ln(x> , (1) 

with = /cbT. A principal virtue of this formulation^ 
is that is subsumes a van der Waals picture of the sol- 
vation without requiring that intermolecular interactions 
of different types be expressed is some specific format. 
We will use this formulation here to distinguish pack- 
ing contributions from longer-ranged interactions that 
are attractive on balance. The notation (...) indicates 
the usual average over the thermal motion of the sys- 
tem. The notation ((. . .))o indicates the average over 
the thermal motion of the system together with an addi- 
tional molecule with no interaction between them, thus 
the subscript 0. ((x))o then the probability that the 
defined inner-shell is empty in the case that interactions 
between the distinguished molecule and the solution are 
absent. That contribution thus gauges the free energy 
cost of finding space for positioning the additional PC 
molecule in the liquid. 

More broadly, these formalities follow from the 
identit>^2 



(x) 



(x) 



(2) 



whic h is a disguised expression of the rule of aver- 
agesW^ This broader observation Eq. might per- 
mit design of weight functions x foi" improved numerical 
performance.^^ 

P{s\x = 1) that appears in Eq. ([T]) is the probabil- 
ity distribution of the binding energies, and is condi- 
tional on an empty inner-shell. That conditioning can 
make P{£\x = 1) simple enough that a Gaussian (or nor- 
mal) approximation suffices .^^^^ The observed behavior 
of P{e) (FIG. fl2| already suggests that possibility. Ac- 
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FIG. 12. Binding energies, centered (be — s — (s)) and 
scaled (cr^ = (^e^)), are approximately normally distributed. 
The slight super-gaussian behavior on the crucial right-side 
of these plots is the signature of repulsive intermolecular in- 
teractions. 



cepting that Gaussian approximation 

/?M(-)p.-ln((x))o 

+ I3{e\x = l) + I3^ {5s'\x = l)/2 

+ In {x> 



(3) 



For conceptual clarity, let us discuss how definition of x 
might be approached. x=l corresponds to the absence of 
van der Waals contact of the distinguished PC molecule 
with the solution. We might chose to define x by assign- 
ing van der Waals radii for all atom types. The results of 
FIGS. 9 and 10, characterizing close atom- atom pair dis- 
tances, and particularly distance ordered contributions 
of FIG. 10 would be directly relevant for that. We would 
adjust those radii assignments to realistically large val- 
ues, while targeting values of (x) not too different from 
one. 

With this conceptual background, Eq. ([3| becomes 
pfi<^^^^ « - In {{x))o + /?{£>+ (fe') /2 . (4) 
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This result is interesting for several reasons. Though it 
is suggestive of a van der Waals treatment, it is not lim- 
ited to a first-order mean-field contribution, and addi- 
tionally the assessment of attractive interaction is made 
on the basis of observations for the case that those inter- 
actions are actually operating. Nevertheless, the packing 
contribution — ln((x))o is typically difficult to calculate 
directly. 

To quantitatively characterize the effort to obtain that 
packing contribution, here we focus on getting an esti- 
mate on the basis of the present methods. With observa- 
tion of the density of the vapor phase (FIG. 13), we can 
get the desired free energy on the saturation curve^^ 



VPvap/ 



(5) 



Our intention (tested below) is to assume that the va- 
por phase is approximately ideal, Pji^^^ ~ 0. The result 
Eq. ([5| applies to the coexisting liquid. For the thermo- 
dynamic states of TABLE [llj at the same temperature 
but slightly different pressures, we apply the correction 



Ap-ln(^) (6) 

.Pvap/ \dpJT VPliq/ 



Of course (3 {dfi/dp)j. Ap = pAp/pn^] at T=300K (TA- 
BLE [ll| this term is about 0.003 and we neglect it. The 
rightmost term of Eq. (|6| extracts the ideal contribution 
to the chemical potential change. Collecting all, we find 
that 



ln((x))o 



In 



Pvap 



/?(e)-/?2(fc2)/2,(7) 



characterizes the net effect of intermolecular excluded 
volumes when the vapor pressure is low. 

Except for the highest temperature considered, the 
ideal estimate of the vapor pressure is low and close to 
the experimental vapor pressure (TABLE [l|), so the as- 
sumption of ideality of the vapor is accurate then. At 
T=600K, the vapor pressure is substantial but the ideal 
estimate of the vapor pressure is still within about 20% 
of the experimental vapor pressure. 

We therefore estimate the packing contribution 
-k^T\n{{x))o ^ 15 kcal/mol at 300K, and 14 kcal/mol 
at 400K. Thus volume exclusion effects contribute to the 
solvation free energies at the level of about 2/3rds of the 
magnitude of the attractive interactions (TABLE [ll]), of 
course with opposite sign, when the vapor pressure is low. 
At 300K, ((x))o ^ 1 X 10-1^ 

Attractive interactions stabilize the liquid of course, 
van der Waals attractions make a larger contribution to 
the mean binding energies (TABLE |llj) than do electro- 
static interactions. Contrariwise, electrostatic contribu- 
tions dominate van der Waals attractions in the variances 
of binding energies. 

If electrostatic interactions are considered solely, then 
at the lowest temperature the mean and variance electro- 
static contributions are roughly in the 2:1 proportion that 




FIG. 13. The dotted results were obtained by the WHAM 
method in order to estimate of the vapor density for the 
T=300K case where that value is not satisfactorily obtained 
from observation of the physical simulation. 



is a symptom of satisfactory Gaussian models of solva- 
tion. If the distributions of binding energies are precisely 
Gaussian, then 



{{e)), = {e)+P{5e^) . (8) 

This follows from the general requirement thatPI' 

P(£) =e-(^-^'^)p(°)(£) (9) 

where P^^^ {e) is the distribution of binding energies for 
the uncoupled case associated with the ((...))o averag- 
ing. Setting ((e))^ = for electrostatic interactions, then 
the 2:1 proportion of mean and variance contributions 
is clear in view of the 1/2 in Eq. Q. At the highest 
temperature considered, the mean and variance contribu- 
tions are roughly equal, so Gaussian models of solvation 
are not supported then. The accuracy of the inference 
of the packing contribution at T=600K is therefore less 
convincing also. 

III. CONCLUSIONS 

Propylene carbonate (PC) does not spread completely 
on graphite, but at ambient conditions it wets with a 
contact angle of 31°. Molecular dynamics simulations 
agree with this contact angle after 40% reduction of the 
strength of graphite-C atom Lennard- Jones interactions 
with the solvent, relative to the models used initially.!^ 
The simulation of a nano-scale PC droplet on graphite 
displays a pronounced layering tendency and an Aztex 
pyramid structure for the droplet. Extrapolation of the 
computed tensions of PC liquid- vapor interface estimates 
the critical temperature of PC accurately to about 3%. 
PC molecules lie ffat on the PC liquid- vapor surface. 
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TABLE 11. Bulk calculations at p = 1 atm, N = 600. Energies are in kcal/mol. kt ^ P (^^^) / (^)- Experimental valueJ^^ 
for A^T at T = 25C are in the range 0.5-0.6 GPa~^. Where not indicated explicitly, estimated statistical uncertainties are less 
than one (1) in the least significant digit given. The thermal expansion coefficient implied by these results is 0.94xlO~^/K 
(experimental valueP 0.845 x 10~^ /K) . 



T(K) 


p (nm 3) 


is) 




KT (GPa-^) 


e/eo 


270 


7.5 


-33.2 


-24.4 


0.3 


74.6 ± 2.6 


300 


7.3 


-32.0 


-24.3 


0.3 


63.4 ± 2.2 


320 


7.2 


-31.2 


-23.6 


0.4 


63.0 zb 2.4 


400 


6.7 


-28.4 


-21.1 


0.5 


55.7 ± 2.2 


600 


5.5 


-22.0 


-15.3 


1.9 


39.4 ± 1.4 



and tend to project the propyl carbon toward the va- 
por phase. Close PC neighbors stack carbonate planes 
with the outer (01) oxygen of one molecule snuggled into 
the positively charged propyl end of another molecule so 
that neighboring molecule dipole moments are approx- 
imately antiparallel. The calculated thermal expansion 
coefficient and the dielectric constants for liquid PC agree 
well with experiment. The distribution of PC molecule 
binding energies is closely Gaussian. Evaluation of the 
density of the coexisting vapor then permits estimation 
of the excluded volume contribution to the PC chemical 
potential, and that contribution is about 2/3rds of the 
magnitude of the contributions due to attractive interac- 
tions, with opposite sign. 
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IV. APPENDIX: METHODS 



aged for 10 ns before a 10 ns production equilibrium tra- 
jectory at 300 K. At 400 K and 600 K, 40 ns equilibrium 
trajectories were obtained. The interfacial tensions were 
assessed (FIG. [6| by differencing interfacial stresses in 
the standard way, averaging through the production tra- 
jectories. 

Configurations were sampled from each trajectory at 
every 0.005 ns for further analysis. The electrostatic con- 
tributions to the binding energies (TABLE [nl) were eval- 
uated with the GROMACS reaction field method. Those 
electrostatic contributions were checked against standard 
Ewald evaluations of electrostatic energies, and the differ- 
ences were typically about 0.1 kcal/mol. The differences 
of generalized reaction field alternatives from the other 
methods were substantially larger, i.e. 1-2 kcal/mol. 

Experimental applications almost always involve a 
racemic mixture of PC. The calculations here were for 
the pure liquid of the enantiomer FIG. [T] This was be- 
cause extensive initial calculations made that choice,!^ 
and no results here are expected to be sensitive to that 
distinction. Explicit checking of a few cases in TABLE [T| 
and TABLE [n] confirm that no results are changed signif- 
icantly for the racemic mixture. Nevertheless, we expect 
subsequent results to treat the racemic case. 



The GROMACS packag^ was used in all simula- 
tions. All the simulations were performed under periodic 
boundary conditions, electrostatic interactions calculated 
by the particle mesh Ewald method with a grid spacing 
of 0.1 nm. Partial charges of PC are those of Ref. 3. cut- 
off at 0.9 nm, and the temperature was maintained by 
the Nose- Hoover thermostat. 

For bulk PC systems, constant NPT conditions were 
adopted, treating 600 PC molecules initially positioned 
uniformly in a (4.4 nm)^ cubic cell utilizing Packmol.^^ 
The initial configuration was energy- minimized, then 
simulations were carried out at T = 270 K, ... 600 K 
with p=l atm, with 1 fs integration time step, and were 
extended to 50 ns with first 10 ns discarded as aging. 

Interfacial characteristics of liquid PC were investi- 
gated by molecular dynamics of two-phase (liquid- vapor) 
systems in slab geometry (FIG. [5|. 600 PC molecules 
were positioned in a 5.3 x 5.3 x 10.3 nm^ cell, T=300, 
400, and 600K. After minimization, these systems were 



A. PC/graphite simulations 

The PC/graphite simulations used 600 PC molecules 
and three layers of graphite. Each layer has 9122 atoms, 
including 266 capping hydrogen atoms at the edges, to 
make a square surface of about 15.5 nm by 15.5 nm. 
A liquid phase configuration of the PC molecules was 
placed near the graphite surface and equilibrated for 5 
ns. All simulations used constant TVN conditions, at 
T=300 K and a cubic (15.5 nm)^ volume. To test how the 
PC/graphite interactions influence interfacial properties, 
a series of simulations with different Lennard- Jones Slj 
values for the PC-graphite interactions were obtained, by 
scaling slj by factors (0.25, 0.4, 0.5, 0.55, 0.6, 0.75 and 
1). Interactions between all PC atoms and the graphite 
carbon atoms were all scaled. Each system was simu- 
lated for 5 ns. For the optimal slj scaling, simulations 
were carried out for 10 ns. Simulations were carried out 
with and without constraints on the graphite surface uti- 
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lizing a harmonic restraint with force constant of 1000 
kJ/(mol nm^). The simulations uses a cut-off value for 
equal to 0.9 nm. Additional simulations were carried-out 
with a cut-off value of 1.2 nm. The PC/graphite con- 
tact angle was found to insensitive to the cut-off value or 
the use of restraints on the surface. The indicated 40% 
reduction in the strength of graphite-C Lennard- Jones 
interactions may be the simplest adjustment that brings 
the simulation results into consistency with the contact 
angle observation. 



were treated. Trajectories of 5 ns/window (after 1 ns of 
minimization and 5 ns of aging) were used to reconstruct 
density profiles. For the 400K calculation, 41 windows of 
0.1 nm uniform spacing were used. 

The harmonic potential energy U = k{z — Zq)'^ win- 
dowing function was employed, with z is the instanta- 
neous distance of the center of mass of the pulled PC 
molecule from center of mass of the slab and zq, the 
designated minimum of [/, identifying the window po- 
sition, k = 4000 kJ/mol/nm^ (for 300K) and A:=6000 
kJ/mol/nm^ (for 400 K). 



B. Droplet silhouette 

The droplet silhouette (FIG. [3| was obtained in the 
following way:^^ On the basis of the simulation data, the 
mass density was binned. A two dimensional grid was 
used in cylindrical coordinates with z perpendicular to 
the graphite plane and r the radial coordinate. The z 
axis contained the centroid of the droplet. The function 



p(r, z) 



Pliq 

2 



1 — tanh 



^o(^) 



d{z) 



(10) 



which acknowledges that the pvap ^ 0, was fit to the 
binned mass density. Thus the width of the interfacial 
profile, d{z), and the position of the interface, ro{z) were 
obtained for each z layer. Since the density exhibits 
distinct layering parallel to z (FIG. [4|, this procedure 
is particularly effective where analysis of the density in 
constant-z slices is natural, for < z < 1.9 nm. 

For z > 1.9 nm, where the < r < 2 nm portion of 
the fluid interface is roughly parallel to the z=constant 
(graphite) surface, constant-r slices of the density were 
similarly fit to the function 



1 — tanh 



d{r) 



(11) 



In view of FIG. [4j this describes the outer (largest r) 
behavior, where the density is decreasing through piiq/2. 
The two approaches give the same values near a crossing 
point r ^2 nm. Combing these two fits gives FIG. [3| 



D. Liquid PC dielectric constant 

The dielectric constant of uniform PC liquid is evalu- 
ated following standard simulation methods.'^^ The cal- 
culation treated 600 PC molecules under periodic bound- 
ary conditions. Resulting values were averaged from 
the 40 ns production trajectories at constant pressure 
of 1 atm. 



C. Windowed sampling for calculations of coexisting 
vapor densities 

The coexisting vapor densities were evaluated by strat- 
ification on the basis of the Weighted Histogram Analy- 
sis Method (WHAM)^^ to concentrate sampling on the 
low densities of the vapor phase. Windowed calculations 
were performed with 600 PC molecules at temperatures 
300K and 400K adopting the methods described above. 
Initial configuration for each window was obtained by 
pulling a PC molecule from the center of mass of the slab 
to deep into vapor phase at constant rate of 0.01 nm/ps. 
In the 300K case, 15 windows, spaced by 0.1 nm in bulk 
liquid phase and 0.05 nm near interface and vapor phase. 
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